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Abstract 

This paper deals with the chaotic oscillator synchronization. A new approach to detect the 
synchronized behaviour of chaotic oscillators has been proposed. This approach is based on the 
analysis of different time scales in the time series generated by the coupled chaotic oscillators. It 
has been shown that complete synchronization, phase synchronization, lag synchronization and 
generalized synchronization are the particular cases of the synchronized behavior called as "time- 
scale synchronization". The quantitative measure of chaotic oscillator synchronous behavior has 
been proposed. This approach has been applied for the coupled Rossler systems. 

PACS numbers: 05.45.Xt, 05.45.Tp 



* Electronic address: laeh@cas.ssu.runnet.rtJ 



1 



Introduction 



Synchronization of chaotic oscillators is one of the fundamental phenomena of nonlinear 



dynamics. It takes place in many physical [1| 



and biological 



BBS 



processes. 



It seems to play an important role in the ability of biological oscillators, such as neurons, to 



act cooperatively 



11 



121. 



There are several different types of synchronization of coupled chaotic oscillators which 
have been described theoretically and observed experimentally 14,Q;I3- The complete 
synchronization (CS) implies coincidence of states of coupled oscillators xi(t) = X2(t), the 
difference between state vectors of coupled systems converges to zero in the limit t ^ oo, |l7l 

inn 

appears when interacting systems are identical. If the parameters of coupled 



chaotic oscillators slightly mismathch, the state vectors are close |xi(t) — X2(t)| ~ 0, but 
differ from each other. Another type of synchronized behavior of coupled chaotic oscillators 
wihth slightly mismatched parameters is the lag synchronization (LS), when shifted in time, 
the state vectors coincide with each other, xi(t + r) = X2(t). When the coupling between 
ator increases the time lag r decreases and the synchronization regime tends to be CS 



oscil 
one 
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23|. The generalized synchronization (GS) 



H 25, Q i 



introduced for drive- 



responce systems, means that there is some functional relation between coupled chaotic 
oscillators, i.e. X2(t) = F[xi(t)]. 

Finally, it is necessary to mention the phase synchronization (PS) regime. To describe 
the phase synchronization the instantaneous phase 0(t) of a chaotic continuous time series is 

16 , 2^, I3| • The phase synchronization means the entrainment 



usually introduced 
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of phases of chaotic signals, whereas their amplitudes remain chaotic and uncorrelated. 

All synchronization types mentioned above are associated with each other (see for de- 
tail 0,0; Q])' t>ut the relationship between them is not completely clarified yet. For each 
type of synchronization there are their own ways to detect the synchronized behavior of 
coupled chaotic oscillators. The complete synchronization can be displayed by means of 
comparison of system state vectors xi(t) and x^ ft), whereas the lag synchronization can be 
determined by means of the similarity function 2l|. The case of the generalized synchro- 
nization is more intricate because the functional relation F[-] can be very complicated, but 
there are several methods to detect the synchronized behavior of coupled chaotic oscillators, 
such as the auxiliary system approach 29^ or the method of nearest neighbors • 
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Finally, the phase synchronization of two coupled chaotic oscillators occurs if the differ- 
ence between the instantaneous phases of chaotic signals ^i,2it) is bounded by some 
constant 

\Mt) - Mt)\ < const. (1) 
It is possible to define a mean frequency of chaotic signal 

n = hm ^ = (0(t)), (2) 

which should be the same for both coupled chaotic systems, i.e., the phase locking leads 
to the frequency entrainment. It is important to notice, to obtain correct results the mean 
frequency Cl of chaotic signal x(t) should coincide with the main frequency Qq = 27r/o of the 
Fourier spectrum (for detail, see j^). There is no general way to introduce the phase for 
chaotic time series. There are several approaches which allow to define the phase for "good" 
systems with simple topology of chaotic attractor (so-called "phase coherent attractor"), 
the Fourier spectrum of which contains the single main frequency /q. 

First of all, the plane in the system phase space may exist, on which the projection of the 
chaotic attractor looks like circular band. For such plane the coordinates x and y can be 
introduced, the origin of this coordinate system should be placed somewhere near the center 
of the chaotic attractor projection. In this case the phase can be introduced as an angle 
in considered coordinate system 0, [s^, but for that all trajectories of chaotic attractor 
projection on the (x, ?/)-plane should revolve around the origin. Sometimes, a coordinate 



13, bal. One can also use the 



transformation can be used to obtain a proper projection 
velocities x and y, if the projections of chaotic trajectories on the plane (x, y) always rotate 
around the origin and in some cases this approach is more suitable . Another way 



to define the phase 0(t) of chaotic time series x{t) is the constructing of the analytical 
sig nal LJ, ^ using the Hilbert transform. Moreover, the Poincare secant surface can be 

nn 

used for the introducing of the instantaneous phase of chaotic dynamical system jU, |27( . 
Finally, the phase of chaotic time series can be introduced by means of the continuous 
wavelet transform js^, but the appropriate wavelet function and its parameters should be 
chosen j^. 

All these approaches give correct results for "good" systems with well-defined phase, but 
fail for oscillators with non-revolving trajectories. Such chaotic oscillators are often called as 
"systems with ill-defined phase" . The phase introducing based on the approaches mentioned 



above for the system with ill-defined phase leads usually to incorrect results f31|. Therefore, 
the phase synchronization of such systems can be usually detected by means of the indirect 



indications 
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37| | and measurements jsj]. 
In this paper we propose a new approach to detect the synchronization between two cou- 
pled chaotic oscillators. The main idea of this approach consists in the analysis of the system 
behavior on different time scales that allows us to consider different cases of synchronization 
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from the universal positions j38|]. Using the continuous wavelet transform 
introduce the continuous set of time scales s and associated with them instantaneous phases 
(f)s{t)- (In other words, is continuous function of time t and time scale s). As we 

will show further, if two chaotic oscillators demonstrate any type of synchronized behavior 
mentioned above (CS, LS, PS or GS), in the time series xi 2(t) generated by these systems 
there are time scales s necessarily correlated for which the phase locking condition 

-0s2(t)| < const (3) 

is satisfied. In other words, CS, LS, PS and GS are the particular cases of the synchronous 
coupled chaotic oscillator behavior called as "time-scale synchronization" (TSS). 

The structure of this paper is the following. In section H] we discuss the continuous 
wavelet transform and the method of the time scales s and associated with them phases 
(f)s{t) definition. In section E] we consider the case of the phase synchronization of two 
mutually coupled Rossler systems. We demonstrate the application of our method and 
discuss the relationship between our and traditional approaches. Section ITTTl deals with the 
synchronization of two mutually coupled Rossler systems with funnel attractors. In this 
case the traditional methods for phase introducing fail and there is no possibility to detect 
the phase synchronization regime, respectively. The synchronization between systems can 
be revealed here only by means of the indirect measurements (see for detail js^). We 
demonstrate the efficiency of our method for such cases and discuss the correlation between 
PS, LS and CS. In section HVl we consider application of our method for the unidirectional 
coupled Rossler systems when the generalized synchronization is observed. The quantitative 
measure of synchronization is described in section El The final conclusion is presented in 
section IVIl 



4 



I. CONTINUOUS WAVELET TRANSFORM 



The continuous wavelet transform is the powerful tool for the analysis of nonlinear dy- 
namical system behavior. In particular, the continuous wavelet analysis has been used for 



the detection of synchronization of chaotic oscillations in the brain 35|, |43|, and chaotic 



laser array 45|. It has also been used to detect the main frequency of the oscillations in 



nephron autoregulation |46[. We propose to analyze the dynamics of coupled chaotic oscil- 
lators using the consideration of system behavior on different time scales s each of them is 
characterized by means of its own phase 0s (t), respectively. So, in order to define the con- 
tinuous set of instantaneous phases (f)s{t) the continuous wavelet transform is the convenient 
mathematical tool. 

Let us consider continuous wavelet transform of chaotic time series x{t) 

+00 

Wis,to)= J xit)rs,toit)dt, (4) 

—00 

where ipsfyit) is the wavelet-function related to the mother-wavelet V'o(^) as 

^..oW = ^^(^)- (5) 



The time scale s corresponds to the width of the wavelet function ips,to{'t)y ^ind to is shift 
of wavelet along the time axis, the symbol "*" in (j3]) denotes complex conjugation. It 
should be noted that the time scale s is usually used instead of the frequency / of Fourier 
transformation and can be considered as the quantity inversed to it. 

r 

The Morlet-wavelet I4j 



1 /-ON f ~V 



^o(^) = -^expijQov) exp (^-^ 



,2 



(6) 



has been used as a mother-wavelet function. The choice of parameter value f^o = 27r provides 
the relation s = 1/ f between the time scale s of wavelet transform and frequency / of Fourier 
transformation. 

The wavelet surface 

W{s,to) = \W{s,to)\e^'^^^'''^ (7) 

describes the system's dynamics on every time scale s at the moment of time to- The value 
of to)| indicates the presence and intensity of the time scale s mode in the time series 



x{t) at the moment of time to- It is possible to consider the quantities 

E{s,to) = \W{s,to)\' 



(8) 



and 

{E{s)) = j \W{sM)?dto, (9) 

which are instantaneous and integral energy distributions on time scales, respectively. 

At the same time, the phase 0s(t) = arg W{s,t) is naturally introduced for every time 
scale s. It means that it is possible to describe the behavior of each time scale s by means 
of its own phase 0s (t). If two interacting chaotic oscillators are synchronized it means that 
in time series xi(t) and X2(t) there are scales s correlated with each other. To detect such 
correlation one can examine the condition Q which should be satisfied for synchronized 
time scales. 



II. PHASE SYNCHRONIZATION OF TWO ROSSLER SYSTEMS 



Let us start our consideration with two mutually coupled Rossler systems with slightly 



mismatched parameters 
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ii,2 = -^i,2yi,2 - ^1,2 + ^(a;2,i - a;i,2), 

2/1,2 = ^i,2a;i,2 + ^2/1,2, (10) 
2:1,2 = P + 2:1,2 (2:1,2 - c), 

where a = 0.165, p = 0.2, and c = 10. The parameters uji^2 = ujq ± A determine the 
parameter mistuning, e is the coupling parameter {uq = 0.97, A = 0.02). In |21i] it has been 
shown that for these control parameter values and coupling parameter e = 0.05 the phase 
synchronization is observed. 

For this case the phase of chaotic signal can be easily introduced in one of the ways 
mentioned above, because the phase coherent attractor with rather simple topological prop- 
erties is realized in the system phase space. The attractor projection on the (x, ?/)-plane 
resembles the smeared limit cycle where the phase point always rotates around the origin 
(Fig. ^a). The Fourier spectrum S{f) contains the main frequency peak /o ^ 0.163 (see 
Fig. ^b) which coincides with the mean frequency / = Q/2tt, determined from the in- 
stantaneous phase (f){t) dynamics Therefore, there are no complications to detect the 
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phase synchronization regime in two coupled Rossler systems fjlOj) by means of traditional 
approaches. 

When the coupling parameter e is equal to 0.05 the phase synchronization between chaotic 
oscillators is observed. The phase locking condition ^ is satisfied and the mean frequencies 
Hi 2 are entrained. So, the time scales Sq — 6 of both chaotic systems corresponding to the 
mean frequencies fii 2 should be correlated with each other. Correspondingly, the phases 
4'si,2{t) associated with these time scales s should be locked and the condition Q should be 
satisfied. The time scales which are the nearest to the time scale Sq should also be correlated, 
but the interval of the correlated time scales depends upon the coupling strength. At the 
same time, there should be time scales which remain uncorrelated. These uncorrelated time 
scales cause the difference between chaotic oscillations of coupled systems. 

Figure |21 illustrates the behavior of different time scales for two coupled Rossler sys- 
tems (fTUj) with phase coherent attractors. It is clear, that the phase difference (psiif) — (ps2it) 
for scales sq = 6 is bounded and therefore time scales sq = 6 corresponding to the main 
frequency of Fourier spectrum /o are synchronized. It is important to note that wavelet 
power spectra {Ei^2{s)) close to each other (see Fig. |2la) and time scales s characterized by 
the large value of energy (e.g., s=5) close to the main time scale Sq = 6.0 are correlated, 
too. There are also time scales which are not synchronized, like s = 3.0, s = 4.0, etc. (see 
Fig. 116). 

So, the phase synchronization of two mutually coupled chaotic oscillators with phase 
coherent attractors manifests itself as synchronous behavior of the time scales sq (and time 
scales s close to sq) corresponding to the chaotic signal mean frequency fl. 

III. SYNCHRONIZATION OF TWO ROSSLER SYSTEMS WITH FUNNEL AT- 
TRACTORS 

Let us consider the more complicated example when it is impossible to correctly introduce 
the instantaneous phase 0(t) of chaotic signal x(t). It is clear, that for such cases the 
traditional methods of the phase synchronization detecting fail and it is necessary to use 
the other techniques, e.g., indirect measurements On the contrary, our approach gives 
correct results and allows to detect the synchronization between chaotic oscillators as easily 
as before. 
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To illustrate it we consider two non-identical coupled Rossler systems with funnel attrac- 
tors (Fig. El): 

il,2 = -^^1,22/1,2 - 2:1,2 + e{x2,l - Xi,2), 

2/1,2 = t^i,2a;i,2 + a?/i,2 + £(2/2,1 - 2/1,2), (H) 
ii,2 = P + 2;i,2(xi,2 - c), 
where e is a coupling parameter, ui = 0.98, uj2 = 1.03. The control parameter values 
have been selected by analogy with js^ as a = 0.22, p = 0.1, c = 8.5. It is necessary 
to note that under these control parameter values none of the methods mentioned above 
permits to define phase of chaotic signal correctly in whole range of coupling parameter e 
variation. Therefore, nobody can determine by means of direct measurements whether the 
synchronization regime takes place for several values of parameter e. On the contrary, our 
approach allows to detect synchronization between considered coupled oscillators easily for 
all values of coupling parameter. 

In it has been shown by means of the indirect measurements that for the coupling 
parameter value e = 0.05 the synchronization of two mutually coupled Rossler systems flllll 
takes place. Our approach based on the analysis of the dynamics of different time scales s 
gives analogous results. So, the behavior of the phase difference (psiit) — 4>s2(t) for this case 
has been presented in figure |3J 6. One can see that the phase locking takes place for the time 
scales s = 5.25 which are characterized by the largest energy value in the wavelet power 
spectra {E{s)) (see Fig. Hla). 

It is important to note that the phase difference (f)si{t) — <ps2{t) is also bounded on the 
time scales close to s = 5.25. So, we can say that the time scales s = 5.25 (and close to 
them) of two oscillators are synchronized with each other. At the same time the other time 
scales (e.g., s = 4.5, 6.0 et. al.) remain uncorrelated. For such time scales the phase locking 
has not been observed (see Fig. HI 6). 

It is clear, that the mechanism of the synchronization of coupled chaotic oscillators is 
the same in both cases considered in sections |n] and IIIII The synchronization phenomenon 
is caused by the existence of time scales s in system dynamics correlated with each other. 
Therefore, there is no reason to divide the considered synchronization examples into different 
types. 

It has been shown in [2l| that there is certain relationship between PS, LS and CS 
for chaotic oscillators with slightly mismatched parameters. With the increase of coupling 
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strength the systems undergo the transition from unsynchronized chaotic oscillations to 
the phase synchronization. With a further increase of coupling the lag synchronization is 
observed. Further increasing of the coupling parameter leads to the decreasing of the time 
lag and both systems tend to have the complete synchronization regime. 

Let us consider the dynamics of different time scales s of two nonidentical mutually 
coupled Rossler systems when the coupling parameter value increases. If there is no 
phase synchronization between the oscillators, then their dynamics remain uncorrelated for 
all time scales s. Figure El illustrates the dynamics of two coupled Rossler systems when 
the coupling parameter e is small enough {e = 0.025). The power spectra {E{s)) of wavelet 
transform for Rossler systems differ from each other (Fig. El a), but the maximum values of 
the energy correspond approximately to the same time scale s in both systems. It is clear, 
that the phase difference (psi(t)—(t>s2(t) is not bounded for almost all time scales (see Fig. 06). 
One can see that the phase difference ipsiit) — ips2(t) increases for time scale s = 3.0, but 
decreases for s = 4.5. It means that there should be the time scale 3 < s* < 4.5 the phase 
difference on which remains bounded. This time scale s* plays a role of a point separating 
the time scale areas with the phase difference increasing and decreasing, respectively. In 
this case the measure of time scales on which the phase difference remains bounded is zero 
and we can not say about the synchronous behavior of coupled chaotic oscillators (see also 
section El . 

As soon as any of the time scales of the first chaotic oscillator becomes correlated with 
the other one of the second oscillator (e.g., when the coupling parameter increases), the 
phase synchronization occurs (see Fig. The time scales s characterized by the largest 
value of energy in wavelet spectrum {E{s)) are more likely to become correlated first. The 
other time scales remain uncorrelated as before. The phase synchronization between chaotic 
oscillators leads to the phase locking Q on the correlated time scales s. 

When the parameter of coupling between chaotic oscillators increases, more and more 
time scales become correlated and one can say that the degree of the synchronization grows. 
So, with the further increasing of the coupling parameter value (e.g., e = 0.07) in the coupled 
Rossler systems the time scales which were uncorrelated before become synchronized 
(see Fig.ini^). It is evident, that the time scales s = 4.5 are synchronized in comparison with 
the previous case {e = 0.05, Fig. El 6) when these time scales were uncorrelated. The number 
of time scales s demonstrating the phase locking increases, but there are nonsynchronized 
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time scales as before (e.g., the time scales s = 3.0 and s = 6.0 remain still nonsynchronized) . 

Arising of the lag synchronization ^] between oscillators means that all time scales 
are correlated. Indeed, from the condition of the lag-synchronization xi{t — r) ~ X2{t) one 
can obtain that Wi{s,t — r) ~ W2{t, s) and therefore (f)si{t — r) ~ <f)s2{t)- It is clear, in this 
case the phase locking condition (jH)) is satisfied for all time scales s. For instance, when 
the coupling parameter of chaotic oscillators (fTT|) becomes large enough (s = 0.25) the lag 
synchronization of two coupled oscillators appears. In this case the power spectra of wavelet 
transform coincide with each other (see Fig. [71a) and the phase locking takes place for all 
time scales s (Fig. 06). It is important to note that the phase difference (psi(t) — (ps2(t) 
is not equal to zero for the case of the lag synchronization. It is clear that this difference 
depends on the time lag r. 

Further increasing of the coupling parameter leads to the decreasing of the time lag r 21 1. 



Both systems tend to have the complete synchronization regime Xi{t) ~ X2{t), so the phase 
difference 0si(t) — 4's2{t) tends to be a zero for all time scales. 

The dependence of synchronized time scale range [sm] Sb] on coupling parameter has been 
shown in Fig. |H1 The range [sm] Sb] of synchronized time scales appears at e ~ 0.039. 
The appearance of synchronized time scale range corresponds to the phase synchronization 
regime. When the coupling parameter value increases the range of synchronized time scales 
expands until all time scales become synchronized. Synchronization of all time scales means 
the presence of lag synchronization regime. 

So, we can say the time scale synchronization (TSS) is the most general synchronization 
type uniting (at least) PS, LS and CS regimes. 

IV. GENERALIZED SYNCHRONIZATION REGIME 

Let us consider another type of synchronized behavior, so-called the generalized syn- 
chronization. It has been shown above, that PS, LS and CS are naturally interrelated 
with each other and the synchronization type depends on the number of synchronized time 
scales. The details of the relations between PS and GS is not at all clear. There are several 
Q) 22 1 dealing with the problem, how GS and PS are correlated with each other. For 



works 



instance, in 



22l | it has been reported that two unidirectional coupled Rossler systems can 



demonstrate the generalized synchronization while the phase synchronization has not been 
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observed. This case allows to be explained easily by means of the time scale analysis. The 
equations of Rossler system are 

xi = -uJiHi - zi, 
iji = uJiXi + ayi, 

zi = P + zi{xi - c) ^^^^ 

X2 = -UJ2y2 - Z2+ e{Xi - X2), 

y2 = UJ2X2 + ay2, 

Z2 =P + Z2{X2 - C), 

where xi = {xi,yi, Zi)'^ and X2 = {x2,y2, 22)"^ are the state vectors of the first (drive) and 
the second (response) Rossler systems, respectively. The control parameter values have 
been chosen as ui = 0.8, uj2 = 1.0, a = 0.15, p = 0.2, c = 10 and e = 0.2. The generalized 
synchronization takes place in this case (see j2^ for detail). Why it is impossible to detect 
the phase synchronization in the system p2|) despite the generalized synchronization is 
observed becomes clear from the time scale analysis. 

Let us consider Fourier spectra of coupled chaotic oscillators (see Fig. IH)). There are two 
main spectral components with frequencies f\ = 0.125 and /2 = 0.154 in these spectra. 
The analysis of behavior of time scales has shown that both the time scales si = I//1 = 8 
of coupled oscillators corresponding to the frequency /i and time scales close to si are 
synchronized while the time scales S2 = I//2 — 6.5 and close to them do not demonstrate 
synchronous behavior fFig. 1101 6). 

The source of such behavior of time scales becomes clear from the wavelet power spectra 
{E{s)) of both systems (see Fig. ^[a). The time scale Si of the drive Rossler system is 
characterized by the large value of energy while the part of energy being fallen on this scale 
of the response system is quite small. Therefore, the drive system dictates its own dynamics 
on the time scale Si to the response system. The opposite situation takes place for the time 
scales S2 (see Fig. [TUla). The drive system can not dictate its dynamics to the response 
system because the part of energy being fallen on this time scale is small in the first Rossler 
system and large enough in the second one. So, time scales S2 are not synchronized. 

Thus, the generalized synchronization of the unidirectional coupled Rossler systems ap- 
pears as the time scale synchronized dynamics, as another synchronization types does before. 
It is also clear, why the phase synchronization has not been observed in this case. Fig. |H1 
shows that the instantaneous phases 0i,2(^) of chaotic signals xi_2(^) introduced by means of 
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traditional approaches are determined by both frequencies /i and /2, but only the spectral 
components with the frequency /i are synchronized. So, the observation of instantaneous 
phases 4>i^2(t) does not allow to detect the phase synchronization in this case although the 
synchronization of time scales takes place. 

Thus, one can see that there is a close relationship between different types of the chaotic 
oscillator synchronization. According to results mentioned above we can say that PS, LS, 
CS and GS are particular cases of TSS. Therefore, it is possible to consider different types 
of synchronized behavior from the universal position. Unfortunately, it is not clear, how one 
can distinguish the phase synchronization. (Here we mean the phase synchronization be- 
tween chaotic oscillators takes place if the instantaneous phase of chaotic signal may be 
correctly introduced by means of traditional approaches and the phase locking condition ^ 
is satisfied), and the generalized synchronization using only the results obtained from the 
analysis of the time scale dynamics. 

V. MEASURE OF SYNCHRONIZATION 

From examples mentioned above one can see that any type of synchronous behavior of 
coupled chaotic oscillators leads to arising of the synchronized time scales. Therefore, the 
measure of synchronization can be introduced. This measure p can be defined as the the 
part of wavelet spectrum energy being fallen on the synchronized time scales 



where [sm] Sb] is the range of time scales for which the condition ((T)) is satisfied and £'1,2 is 
a full energy of the wavelet spectrum 



This measure p is for the nonsynchronized oscillations and 1 for the case of complete 
and lag synchronization regimes. If the phase synchronization regime is observed it takes 
a value between and 1 depending on the part of energy being fallen on the synchronized 
time scales. So, the synchronization measure p allows not only to distinguish the synchro- 
nized and nonsynchronized oscillations, but characterize quantitatively the degree of TSS 
synchronization. 




(13) 




(14) 
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Fig. presents the dependence of the TSS synchronization measure pi for the first 
Rossler oscillator of system (fTT|) considered in section IIIII on the coupling parameter e. It 
is clear that the part of the energy being fallen on the synchronized time scales growths 
monotonically with the growth of the coupling strength. Similar results have been obtained 
for the generalized synchronization of two coupled Rossler systems considered in section IIVI 

It has already mentioned above that when the coupled oscillators do not demonstrate 
synchronous behavior there are time scales s* the phase difference (psi{t) — Vs2{t) on which 
is bounded. Such time scales play role of points separating the time scale areas where the 
phase difference increases and decreases, respectively (see also section ITTT|l . Nevertheless, the 
presence of such time scales does not mean the occurrence of chaotic synchronization because 
the part of energy being fallen on them is equal to zero. Therefore, the synchronization 
measure p of such oscillations is zero, and dynamical regime being realized in the system in 
this case should be classified as non-synchronous. 

VI. CONCLUSION 

Summarizing this work we would like to note several principal aspects. Firstly, we have 
proposed to consider the time scale dynamics of coupled chaotic oscillators. It allows us 
to consider the different types of behavior of coupled oscillators (such as the complete syn- 
chronization, the lag synchronization, the phase synchronization, the generalized synchro- 
nization and the nonsynchronized oscillations) from the universal position. In this case TSS 
is the most common type of synchronous coupled chaotic oscillator behavior. Therefore, 
the another types of synchronous oscillations (PS, LS, CS and GS) may be considered as 
the particular cases of TSS. The quantitative characteristic p of the synchronization mea- 
sure has also been introduced. It is important to note that our method (with insignificant 
modifications) can also be applied to dynamical systems synchronized by the external (e.g., 
harmonic) signal. 

Secondly, the traditional approach for the phase synchronization detecting based on the 
introducing of the instantaneous phase (j){t) of chaotic signal is suitable and correct for such 
time series characterized by the Fourier spectrum with the single main frequency /q. In this 
case the phase 0so associated with the time scale Sq corresponding to the main frequency 
/o of the Fourier spectrum coincides approximately with the instantaneous phase (f){t) of 
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chaotic signal introduced by means of the traditional approaches (see also 3^). Indeed, as 
the other frequencies (the other time scales) do not play a significant role in the Fourier 
spectrum, the phase (j){t) of chaotic signal is close to the phase (psoit) of the main spectral 
frequency /o (and the main time scale Sq, respectively). It is obvious, that in this case the 
mean frequencies / = (0(t))/27r and fso = {(psoit)) /2tt should coincide with each other and 
with the main frequency /o of the Fourier spectrum (see also [31]) 



f — fso — fo- 



(15) 



If the chaotic time series is characterized by the Fourier spectrum without the main single 
frequency (like the spectrum shown in the Fig. El 6) the traditional approaches fail. One has 
to consider the dynamics of the system on all time scales, but it is impossible to do it by 
means of the instantaneous phase 0(t). On the contrary, our approach based on the time 
scale dynamics analysis can be used for both types of chaotic signals. 

Finally, our approach can be easily applied to the experimental data because it does 
not require any a-priori information about the considered dynamical systems. Moreover, in 
several cases the influence of the noise can be reduced by means of the wavelet transform 



(for detail, see 



39, 



We believe that our approach will be usefu 



the analysis of physical, biological, physiological and other data, such as 



and effective for 
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CAPTIONS 

Fig. 1. (a) Phase coherent attractor and {b) spectrum of the first Rossler system ()lUp. 
Coupling parameter e between oscillators is equal to zero 

Fig. 2. (a) Wavelet power spectrum {E(s)) for the first (solid line) and the second 
(dashed line) Rossler systems (fTUI). (6) The dependence of phase difference (psi(t) — 4's2if) 
on time t for different time scales s. The coupling parameter between oscillators is £ = 0.05. 
The phase synchronization for two coupled chaotic oscillators is observed 

Fig. 3. (a) Phase picture and (6) power spectrum of the first Rossler system (|TT|l oscil- 
lations. Coupling parameter e is equal to zero 

Fig. 4. (a) The normalized energy distribution in wavelet spectrum {E{s)) for the first 
(the solid line denoted "1") and the second (the dashed line denoted "2") Rossler sys- 
tems (fTTj): (6) the phase difference (psiif) — (psiif) for two coupled Rossler systems. The 
value of coupling parameter has been selected as £ = 0.05. The time scales s = 5.25 are 
correlated with each other and the synchronization has been observed 

Fig. 5. (a) The normalized energy distribution in wavelet spectrum {E{s)) for the first 
(the solid line denoted "1") and the second (the dashed line denoted "2") Rossler systems; (fe) 
the phase difference (psiit) — (f)s2it) for two coupled Rossler systems. The value of coupling 
parameter has been selected as e = 0.025. There is no phase synchronization between 
systems 

Fig. 6. (a) The normalized energy distribution in wavelet spectrum {E{s)) for the first 
(the sohd line denoted "1") and the second (the dashed line denoted "2") Rossler systems; 
(b) the phase difference 4>si(t)—4>s2(t) for two coupled Rossler systems. The value of coupling 
parameter has been selected as e = 0.07. 

Fig. 7. (a) The normalized energy distribution in wavelet spectrum {E{s)) for the Rossler 
system; (6) the phase difference 4>si(t) — 4>s2(t) for two coupled Rossler systems. The value 
of coupling parameter has been selected as e = 0.25. The lag synchronization has been 
observed, all time scales are synchronized 

Fig. 8. The dependence of the synchronized time scale range [sm] Sb] on the coupling 
strength e for two mutually coupled Rossler systems (fTT|) with funnel attractors 
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Fig. 9. Fourier spectra for(a) tlie first (drive) and (b) the second (response) Rosier 
systems fll2|) . Tlie coupling parameter is e = 0.2. The generalized synchronization takes 
place 

Fig. 10. (a) The normalized energy distribution in wavelet spectrum {E{s)) for the first 
(the solid line denoted "1") and the second (the dashed line denoted "2") Rossler systems. 
The time scales pointed with arrows correspond to the frequencies /i = 0.125 and /2 = 0.154, 
respectively; (6) the phase difference 4>siit) — 4>s2it) for two coupled Rossler systems. The 
generalized synchronization has been observed 

Fig. 11. The dependence of the synchronization measure pi for the first Rossler sys- 
tem (fTTj) on the coupling strength e. The measure p2 for the second Rossler oscillator 
behaves in a similar manner, so it has not been shown in the figure 
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